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We extract the string susceptibility exponent j s tr by measuring the distribution of baby universes on surfaces 
in the case of various matter fields coupled to discrete 2d quantum gravity. For c < 1 the results are in good 
agreement with the KPZ-formula, if logarithmic corrections are taken into account for c = 1. For c > 1 it is not 
as clear how to extract j s tr but universality with respect to c is observed in the fractal structure. 



1. Introduction 

The fractal structure of matter coupled to dis- 
crete 2d gravity has been studied extensively dur- 
ing the last two years [3]. The motivation has 
been twofold. First to observe the crossover to 
a branched polymer phase that is expected when 
large amount of matter is coupled to gravity, sec- 
ondly the hope has been to observe some change 
in the fractal structure at c = 1 which could 
explain the breakdown of the continuum formal- 
ism. Until now the simulations have focused on 
measuring things like the average curvature, the 
branching of surfaces and the maximal distance 
between two vertices. These measurements have 
given some indications of universality in the frac- 
tal structure with respect to c but the results have 
not been convincing. This is though not so alarm- 
ing as there are a priori no reasons why these 
quantities should be universal. At the c = 1 bar- 
rier no evidence of any pronounced change in the 
fractal structure has been observed. 

Here we want to report on direct measurements 
of the string susceptibility 7 using a new method 
where we study the distribution of "baby uni- 
verses" on the surfaces. From the KPZ-formula 
we have analytical predictions for 7 for theories 
with c < 1, but for c > 1 it yields complex and 
unphysical values, indicating a disease in the con- 
tinuum formalism. 

The theories we study are multiple copies 
of spin models coupled to dynamical triangu- 
lated random surfaces [1], more explicit: (/-state 
Potts models (q = 2,3,4, with respective cen- 



tral charges c = 1/2, 4/5, 1) and Gaussian models 
(c = 1). The partition function is 

n m 

ZN T ,nJfl)=J2I[ Z >( T ) « 

T k = l 

where the sum is over triangulations T and Zf(T) 
is the partition function for a single model, i.e. 

Z k P (T)= J2 e^<^> S( ° i(k) '° j(fc)) (2) 




For details of the simulations we refer to [1]. 



2. Baby universes and j s t r 

Lets start by a brief description of the method. 
We measures the fractal structure of the surfaces 
by looking at the distribution of baby universes. 
A baby universe is defined as a simply connected 
region (of size B) of a surface (of size Nt) where 
the boundary length / is much smaller than the 
square of the area it encloses (/ <C y/B). 

We are mainly interested in baby universes 
with minimal necklength, mimbu, (as those are 
easiest to identify on the lattice). On a triangu- 
lated surface this means / = 3. We can calculate 
the distribution of mimbu njy T (B) bye asking in 
how many ways it is possible to glue together two 
random surfaces, of sizes B and Nt — B, along 
the boundaries of one triangle. This gives [4] 

n NT (B)~(N T -By- 2 B-<- 2 (4) 



2 



(For arbitrary necklength there is also an / depen- 
dence). This follows from the asymptotic behav- 
ior of the partition function 



Z 



(5) 



The distribution njy T (B) is easily measured and 
7 is found by fitting it to (4). 

But as eq. (5) is only asymptotically correct 
there will be some deviation from (4) for small 
B. In order to compensate for those "finite size" 
effects we have done two things: (1) we included 
a correction term c\jB and fitted to the form 

In n NT (B) = k + (7 - 2) ln[(7V T - B)B] + Cl /B(6) 

(2) we introduced a lower cutoff Bq in the data 
and, in order to treat all the data consistently, 
assumed that the effects of Bq could be approxi- 
mated as 



7b = 7 - c 2e 



-c 3 B 



(7) 



To test this method we have applied it to the 
case of pure gravity (where 7 = —1/2) and the 
results for different lattice sizes gave 



7 



1000 -0.496 ± 0.005 
2000 -0.501 ± 0.004 
4000 -0.504 ± 0.004 



3. Matter fields coupled to gravity 
3.1. c < 1 

Now we would like to apply this to theories 
including matter fields. Lets start with c < 1. 
Then we have looked at two simple spin models; 
one Ising model (c = 1/2, 7 = —1/3) and one 
q = 3 Potts model (c = 4/5, 7 = —1/5). 

The fitted values of 7 can been seen in fig. 1 . 
We get (as expected) the pure gravity value 7 = 
— 1/2 away from the phase transitions, but in the 
vicinity of (3 C we get sharp peaks. And the peak 
values are in good agreement with 7 predicted 
by the KPZ-formula. That we see sharp peaks 
is also in agreement with the expectations that 
these models only couple weakly to gravity (for 
the Ising model 7 is only changed at (3 C [2]). 
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Figure 1. Fitted values of 7 for one Ising model 
and one q = 3 Potts model coupled to 2d gravity. 
( Ay = 1000.) 



3.2. c = 1 

More interesting are theories with c = 1. In 
this case the KPZ-formula predicts 7 = 0, but 
we also know from analytical calculations that 
we should include logarithmic corrections to the 
asymptotic form of Zjy T [4]. Then the expected 
distribution of baby universes (4) is modified and 
we fit to 



In riN T (B) 



-- k + (7 - 
-a[ln(7Vy 



2)ln[(7V T ■ 
- B) In B] 



B)B] 
Ci/B 



(8) 



where a is an unknown parameter. 

We looked at two models that have c = 1, one 
q = 4 state Potts model and one Gaussian model. 
In fig. 4 we show the results for the former, both 
with and without log. corrections. As before we 
get 7 = —1/2 away from the transition and a 
peak around (3 C . If we include the log. correction 
the peak value agrees well with 7 = 0, but we get 
different value otherwise. 

The same thing was seen for one Gaussian 
model coupled to gravity. For Nt — 4000 and 
with log. corrections we got 7 = —0.09 zb 0.08 
(compared to 7 ~ —0.3 without log. correction). 

3.3. c > 1 

For c > 1 we looked at 2 and 4^ = 3 Potts 
models (c = 1.6 and 3.2) and 2 to 5 Gaussian 
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Figure 2. Fitted values of 7 for one q = 4 state 
Potts model coupled to 2d gravity. (Nt = 1000.) 



Figure 3. Fitted values of 7 vs c for all the models 
studied, both with (filled dots) and without (open 
dots) log. corrections. (Nt = 1000.) 



models (c = 2 to 5). For these models we do 
not know which kind of correction term (if any) 
should be included, and, unfortunately, our expe- 
rience in the c = 1 case has shown us that the 
exact functional form is important in order to ex- 
tract the correct value of 7. 

What we have done is to make fits to the dis- 
tributions njy T (B) with and without log. cor- 
rections (with a as free parameter). The result- 
ing peak values are shown in fig. 3 (as a function 
of c) for all the models we looked at. The ex- 
tracted values of 7 fall on two curves, depending 
on weather log. corrections are included or not. 
But the curves are very different. Although this 
shows us that we cannot predict the value of 7 
without some additions knowledge of what kind 
of corrections to include, it is clear from fig. 3 that 
we can claim universality with respect to c as the 
values fall on the same curve regardless if we look 
at multiple Potts or Gaussian models. 

4. Discussion 

From the results for c < 1 we see that this 
method of measuring 7 works very well and yields 
results that are in good agreement with predicted 
values. But the c = 1 models also show us that 
it is crucial to know the exact functional form 
to which the distributions are fitted in order to 
extract correct value of 7. 



This is unfortunately not known in the case c > 
1 and we get different result for 7 depending on 
what kind of correction we include. So we cannot 
claim to have extracted 7 in this case. What we 
can see from the distributions of baby universes 
is that using the same kind of corrections yields 
the same 7 for different models with the same 
central charge, hence we observe universality in 
the fractal structure with respect to c. 
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